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Abstract 

We consider a general, classical theory of gravity with arbitrary 
matter fields in n dimensions, arising from a diffeomorphism invariant 
Lagrangian, L. We first show that L always can be written in a "man- 
ifestly covariant" form. We then show that the symplectic potential 
current (n — l)-form, 0, and the symplectic current (n — l)-form, u;, 
for the theory always can be globally defined in a covariant manner. 
Associated with any infinitesimal diffeomorphism is a Noether cur- 
rent (n — l)-form, J, and corresponding Noether charge (n — 2)-form, 
Q. We derive a general "decomposition formula" for Q. Using this 
formula for the Noether charge, we prove that the first law of black 
hole mechanics holds for arbitrary perturbations of a stationary black 
hole. (For higher derivative theories, previous arguments had estab- 
lished this law only for stationary perturbations.) Finally, we propose 
a local, geometrical prescription for the entropy, Sd yn , of a dynamical 
black hole. This prescription agrees with the Noether charge formula 
for stationary black holes and their perturbations, and is independent 
of all ambiguities associated with the choices of L, 0, and Q. How- 
ever, the issue of whether this dynamical entropy in general obeys a 
"second law" of black hole mechanics remains open. In an appendix, 
we apply some of our results to theories with a nondynamical metric 
and also briefly develop the theory of stress-energy pseudotensors. 
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1 Introduction 



Recently, many authors have investigated the validity of the first law of black hole mechanics 
and the definition of the entropy of a black hole in a wide class of theories derivable from a 
Hamiltonian or Lagrangian In particular, in Q the first law was proven to hold in 

an arbitrary theory of gravity derived from a diffeomorphism invariant Lagrangian, and the 
quantity playing the role of the entropy of the black hole was identified as the integral over the 
horizon of the Noether charge associated with the horizon Killing vector field. Although some 
key issues concerning the validity of the first law and the definition of black hole entropy in a 
general theory of gravity were thereby resolved, the analysis of Q, nevertheless, was deficient 
in the following ways: (1) It was not recognized that a diffeomorphism covariant choice of 
the symplectic potential current form always can be made. Consequently, several steps in the 
arguments were made in an unnecessarily awkward manner. (2) While a completely general 
proof of the first law of black hole mechanics was given for perturbations to nearby stationary 
black holes, a proof of the first law for non-stationary perturbations was given only for theories 
in which the Noether charge takes a particular, simple form. (3) A proposal was made for 
defining the entropy of a dynamical black hole. However, this proposal made use of a rather 
arbitrary choice of algorithm for defining the symplectic potential current form, and it turns 
out to possess the undesireable feature that the addition of an exact form to the Lagrangian 
(which has no effect upon the equations of motion of the theory) can induce a nontrivial change 
in this proposed formula for the entropy of a dynamical black hole [11]. 

The main purposes of this paper are to remedy all of the above deficiencies, and, in addition, 
develop further the theory of Noether currents and charges in diffeomorphism invariant theories. 
We shall show, first, that the Lagrangian of a diffeomorphism invariant theory always can be 
expressed in a manifestly covariant form. This will enable us to give globally defined, covariant 
definitions of the symplectic potential current form, 0, and symplectic current form, in an 
arbitrary diffeomorphism invariant theory. Furthermore, results on the general form of in 
an arbitrary theory will be obtained, from which it will follow that the Noether charge form, 
Q, always has a particular, simple structure. As a consequence of this structure of Q, the 
first law of black hole mechanics will be proven to hold for nonstationary perturbations in an 
arbitrary theory of gravity derived from a diffeomorphism covariant Lagrangian. We then shall 
propose a definition of the entropy, Sd yn , of an arbitrary cross-section of a nonstationary black 
hole, wherein Sdyn is given by an integral over the horizon of a local, geometrical quantity. Our 
proposed definition agrees with the known answer (as determined by the first law) for stationary 
black holes and their perturbations, and is independent of all ambiguities associated with the 
choices of L, 0, and Q. However, it is not known whether our Sd yn obeys a "second law" 
in general theories of gravity. The paper concludes with an Appendix in which some of our 
results are applied to theories with a nondynamical metric, and some results on stress-energy 
pseudotensors are obtained. 
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We shall follow the notation and conventions of |12| . All spacetimes, tensor fields, and 
surfaces considered in this paper will be assumed to be smooth (C°°). 

2 The form of the Lagrangian for diffeomorphism invariant the- 
ories 

We wish to consider, here, Lagrangian theories on an n-dimensional, oriented manifold M, 
with the dynamical fields consisting of a Lorentz signature metric g a b, and other fields ip. For 
simplicity and definiteness, we shall restrict consideration to the case where ip is a collection of 
tensor fields on M (with arbitrary index structure). However, we foresee no essential difficulty 
in extending our analysis and results to the case where ip is a section of an arbitrary vector 
bundle which possesses a connection uniquely determined by g a b- 



We start with the general form of a Lagrangian postulated in [14] and ||: Specifically, we 

o 

introduce an arbitrary, fixed, globally defined, derivative operator, V, and take the Lagrangian 
to be a function of the quantities g a b, ip, and finitely many of their symmetrised derivatives with 

o 

respect to V- In addition, the Lagrangian is permitted to depend on additional "background 

o o 

fields" 7 - which, like V, do not change under variation of the dynamical fields; a good example 

o 

of such a background field upon which the Lagrangian could depend is the curvature, Rbcd e , of 

o 

V- Thus, we take the Lagrangian to be an n-form locally constructed - in the precise sense 
explained in [O] - out of the following quantities, 



L = L {g ab , Vai 9ab, V( ai ••• V flfc ) 9ab, Ip, Voj 4>, •••> V( ai ••• V a; ) i>, 7 J , (1) 

Here and in what follows, we use boldface letters to denote differential forms on spacetime, and 
we shall, in general, suppress their tensor indices. In the following, we also shall collectively 
refer to the dynamical fields a ip and g v as "</>" . 

We shall be concerned here only with diffeomorphism invariant theories, i.e., the Lagrangian 
will be assumed to be diffeomorphism covariant in the sense that 

l(/*(0)) = rm, (2) 

where /* is the action induced on the fields by a diffeomorphism / : M — > M . Note that on the 

o o 

left side of this equation /* does not act on V or the background fields 7- 

The main result to be established in this section is that any Lagrangian, L, which is diffeo- 
morphism covariant in the sense of eq.(Q) always can be written in a manifestly covariant form. 
More precisely, we have the following lemma, which is closely related to "Thomas replacement 
theorem" [15]: 



Lemma 2.1 7/L as given in fy) is diffeomorphism covariant in the sense of ^) then L can be 
re-expressed as 

L = L (g ab ,V ' ai Rbcde, ■-, V( 01 ...V aro )i4 C( fe ! '!/>, V ai rp, V( ai ...V a! )i/>) (3) 

where V denotes the derivative operator associated with g a b, m = max(/c — 2, 1 — 2), R a bcd denotes 
the curvature of g a b, o,nd the absence of any dependence on "background fields" in (Qj should be 
noted. 
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Proof. We begin by using the relation (written here schematically) 

o o 

V ol = Va + a ■ terms linear in V g 

o 

for any tensor field, a, to re-write all of the V-derivatives of the matter fields, tp, in terms 
of V-derivatives of -0 ~ where V is the derivative operator associated with g - together with 

o _ 

terms involving the V-derivatives of g. Next, we re-write the V-derivatives of tp in terms of 
symmetrized V-derivatives and the curvature of g and its derivatives. Then we re-write the 

o o 

curvature of g and its derivatives in terms of V-derivatives of g and the curvature of V and 

o o o 

its V-derivatives. Finally, we write all of the V-derivatives of g in terms of symmetrized V- 

o o 

derivatives of g and the curvature of V and its V-derivatives. We thereby obtain 

(o o o o \ 

g, Voi gab,-, V(oi - V 0s ) 5a6,V , »V ai ^,-,V( 01 ...V 0| )^,7 J (4) 

O O O 

where s = max(A:, I) and 7 is comprised by 7 together with the curvature of V and (finitely 

o o o 

many of) its V-derivatives. Next we eliminate Va gbc and its higher V-derivatives in favor of 

C e cd = ^g ef {V c gfd+Vdgf C -Vf g C d) (5) 

o 

and its V- derivatives via the substitution 

o 

Va gbc = gecC e ab + gbeC £ ac (6) 
o o 

Again, we express all V-derivatives of C e c a in terms of symmetrized V-derivatives and the 

o 

curvature of V- We thereby obtain 

L = L C e cd , Voi C e cd , V (oi - V ._!) C e cd , 0, V ai V, V (ai ...V ai) V, 7 ') . (7) 

It is tedious but straightforward to check that the symmetrized derivatives of C can be re- written 

as 

o o o o 

V(ai ••• Va p ) C 6 cd = V(ai ••• Va p C^cd) 

p + 3 



+ 4(J3+ l)(p + 2) ^ V ( a l-"^ a '-" Va p) (-^ca^ 6 



3p + 4 



^(V (d V ai ...V ai V a ,...V ap) i? ai 



8p(p + l)(p + 2)^ 
+V (c V ai ...V ai V a ....V ap) i? ajda . e ) 

o 

+terms involving no more than (p — 1) V- derivatives of C (8) 
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where V fli means the omission of this derivative operator in the sequence. By repeatedly making 
this substitution in sequence, starting with p = s — 1, then p = s — 2, etc., and, at each step, 
writing multiple derivatives in terms of symmetrized derivatives and curvatures, we can express 
the Lagrangian as 



L — L ydab, C e c d, V(ai C 6 cd)) ■•■> V(ai ••• Va s _i C & cd) i Rbcdei V ai -R& C( fe> V( 01 ---^ a B _ 2 )Rbcdei 

^,V ai ^,...,V (ai ...V ai) V,7') (9) 



The infinitesimal version of the diffeomorphism covariance condition (0) is 



(9/ 



Applying this to eq.(EJ), we obtain 



<9L <9L o 
-^——L^C cd H 5 L € V( aa ^ cd) + ••• 

dC cd 9 v (oi C e crf) 

dli ° ° <9L o t 

H 5 5 ^5 ••• Va s _! C cd) H — ^ 7 

t^t^oC cd + — V(ai oC cd) + ... 

<9L o o 

+ —o 5 V (oi - Va s _! SC e cd) (11) 



5C° 6c = g ad ( V( 6 V C ) & - ^ (fec) e e e ) - 2 V (a £ d W^ c (12) 



where 



is the variation of C a b c arising from the metric variation 5g a b = C^gab = 2V( ^)- (Note that the 
terms in eq.(|lC)|) arising from the variations of ip, g, and the curvature, R, of g cancel and were 
therefore omitted when writing eq.([Tl|)). The dependence of the terms in eq. (0) on £ a and its 

o 

symmetrized V-derivatives should be noted: no more than one derivative of £ a appears on the 
left side of this equation, but the right side contains terms with as many as (s + 1) symmetrized 
derivatives of £ a . Since, at any given point in M, £ a and its symmetrized derivatives can be 



chosen independently, it follows directly that a necessary condition for eq.(ll) to hold for all £ a 
is 

—5 = 0, for i = 0...(s - 1). (13) 

V (oi .» Va, C e cd) 

This reduces the lagrangian to the form 



L = L (g ab ,R bcde ,V ai Rbcde,---,^( ai ---^a m )Rbcde,4',^a 1 ^,---,^(a 1 ---^a l )^,^j , (14) 



where m = s — 2 = max(/c — 2,1 — 2). The diffeomorphism invar iance condition (|il|) yields one 
more relation, namely, 

-^C C 7' = (15) 

Qry, 
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where a sum over the fields 7 ' should be understood. To show that this implies that L has no 

o 1 

essential dependence on 7 , we proceed by introducing a local coordinate system x , x n , and 

o 

viewing L as a function of the coordinate components of the dynamical fields and 7 . We then 

o 

view the components of 7 as given functions of x^. In this way, we may view L as 



L = L 



(5, Rbcde, ^ ai Rbcde, ■■; V( 0j . . .V am) R bcde , lp , V 0l ^, V( ttl ...V aj )V>, X^j (16) 



i.e., we replace the dependence of L on 7 by explicit dependence on coordinates. Condition 
( p~T| ) then implies that 

V^" = ' ( 17 ) 

Clearly, this equation holds for all £ a if and only if 

dL 



dx^ 



0. (18) 



We therefore see that ( |18D implies that any diffeomorphism invariant lagrangian must be of the 
form 

L = L [g a b, Rbcde, V ai R bcde ,..., V( ai ...V am )i4 C( fe,^, V ai V>, V( ai ...V aj )^J . (19) 
as we desired to show. □ 

3 The form of the symplectic potential and symplectic currents 
for diffeomorphism invariant theories. 



As is well known (and as will be explicitly demonstrated in Lemma |3.1| below), if we vary the 
dynamical fields = (g a ft,'0), then - by "integration by parts" manipulations of the terms 
involving the derivatives of 5(f) - the first variation of the Lagrangian can always be expressed 
in the form 

5L = E5(f> + d@ (20) 

with 

E8<t> = (E g ) ab 5g ab + E^V (21) 

where a sum over the "matter fields" ifi is understood, and it is also understood that for each 
matter field, has tensor indices dual to ip, and these indices are contracted with those of 
Sip in eq.(21). Here E 9 and are locally constructed out of the dynamical fields <f> and their 



derivatives, whereas is locally constructed out of <fi, 5(f) and their derivatives and is linear in 
5(f). The equations of motion of the theory are then taken to be 

(E g ) ab = , and E^ = (22) 

The (n — l)-form, 0, defined by eq.(pO|) is called the symplectic potential form (see below). 
However, although the equations of motion form, E, is uniquely determined by eq.(^), this 
equation determines only up to the addition of a closed (and hence exact p^| ) (n — l)-form. 
Thus, some arbitrariness is present in the choice of 0. The principal result of this section is 
stated in the following lemma: 
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Lemma 3.1 Given a covariant lagrangian of the form (\19{ ) one can always choose a covariant 
satisfying (ffiy. Moreover, can be chosen to have the form 



= 2E b f V d Sg bc + ©' (23) 

where 0' is of the form 

m-l l-l 

©' = S ab (cf>)S 9ah + ^T i (0) a6c ^'-- a ^V (ai ...V Oi) i? a6cd + ^U l <«/»)^--^<5V (ai ...V Oi) ^ (24) 

i=0 i=0 

In other words, in the expression for 0, the 5's can be put to the left of derivatives of the 
dynamical fields everywhere except for the single term E b ^ d V d 5g bc ■ Finally, is given by 

{&R J& 2 ...5n - ^R e ab 2 ...b n [*0) 

where E R bcd e blb2bn is the equation of motion form that would be obtained for R abcd if it were 
viewed as an independent field in the Lagrangian ftlQ ) rather than a quantity determined by the 
metric. 



Proof. Given L in form ( |19| ) we write it as L = Le, where e is the canonical volume form on M 
associated with g ab . Computing a first variation, we obtain 

(8L 8L 8L 8L 
-7T— + T7T> ^ Rabcd + 7^ 5 SV ai R abcd + ... + — (5V (ai ■ ■ -V am) R abcd 
9a b Oft-abed d\/ ai rC a bcd & V (01 ••■ ^ a m ) K abcd 

81 81 81 \ 1 

(For tensors, such as 5R abcd , whose components are not algebraically independent at each point, 
we uniquely fix the partial derivative coefficients appearing in this equation by requiring them 
to have precisely the same tensor symmetries as the varied quantities.) In order to obtain the 
desired expression for 0, we must suitably rewrite eq.(26) in the form (pQ). To see how this 
can be done, we focus attention on a typical term, 

8L 



^V( ai ...V a .)-R a jcd 

and rewrite it as, 



SV {ai ...V ai) R abcd (27) 



8L ( 8L 

SV( ai ...V ai ) Rabcd = 6 ^ p V ai (<5V a2 ...V ai .Rabcd) 



^V( ai ...V a .)i? a 6cd 1 "' ' \5V( ai ...V 0i )i?a6cd 



+e • terms proportional to V5g 

V «i ( e ^v7 V7 n $Va 2 ---Va l Rabcd J 

\ , V( ai ...V ai )ita6cd / 

+V ai (e • (terms proportional to Sg)) 

~ V ai ( e ^ ^ B I ^a 2 ---^a t Rabcd 

\ 0\/( ai ...V a .)K a bed j 

+e • (terms proportional to bg) 

dV - eV ai I— ^— ] 6V a2 ...X7 ai R abcd 

+e ■ ( terms proportional to 5g) (28) 
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where the (n — l)-form, V, has the form 

dL 

V b2 ... K — e aib2...b n tj^ ~ — — <5V a2 ...V ai i? a fe C£ i — e • (terms proportional to Sg) (29) 

OV (a 2 ---V ai ) Rabcd 



This shows that we can rewrite our original term (|27| ) as a sum of a similar term of lower differ- 
ential order, the exact form dV, and terms proportional to 5g. By iterating this procedure and 
performing similar manipulations on all other terms in (|26|) containing derivatives of variations 
of the curvature or matter fields, we obtain, 

5L = e (Af5g ab + E abcd 5R abcd + E^) + d® (30) 



where the (n — l)-form, 0, is covariant and has the same structure as the right side of eq.(24). 
Note that in eq.(^) eE^ are precisely the equations of motion form for the matter fields ip and 

£<C " = 6 (alb " Vai dvBu^ + - + (- 1 ) mv K--- v -) 9V(ai ...v am) ^ c J ( 31 ) 

would be the equations of motion form for R a bcd ^/it were viewed as an independent field. In 
fact, however, R abcd is n °t an independent field, and, taking account of the symmetries of E R bcd , 
we have 

E a R bcd 5R abcd = 2E abcd V a V d 5g bc + E a R bcd R abc e Sg de (32) 
Making this substitution and integrating twice by parts, we obtain 

5L = e(A b g c 5g bc + 2V a V d E R bcd 5g bc + E^) (33) 

+d(2E bcd V d g bc - 2V d B bcd 5g bc + 0) (34) 

where 

fT?bcd\ rpabcd^ /q r \ 

O )b 2 ...b n - E R e ab2 ,„ bn (35) 

and 

Note that the equations of motion associated with g ab are thus 

\bc / Jfec I qt-7 v7 j^a,bcd\ 



A bc = A bc + E™ rb R pqr c . (36) 



(E g ) bc = e + 2V a V d E a R bcd ) ■ (37) 
Thus, we obtain 

= 2E b £ d V d 5g bc + 0', (38) 

where 

0' = - 2V d E% d 5g bc . (39) 

This shows that is manifestly covariant and of the form claimed in the statement of the 
lemma. □ 

We comment, now, on the possible ambiguities in the choice of for a covariant Lagrangian. 
The above lemma proves that always can be chosen to be covariant. This appears to be a 
very natural requirement, and, in the following, we shall restrict consideration to covariant 
choices of 0. The statement of the lemma also provides the canonical form (23) for 0, which 
will play an important role in our analysis below. However, this general form does not uniquely 
determine 0, since one could add to an exact (n — l)-form which has the structure of the 
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the right side of eq.(24). The proof of the lemma does implicitly provide a particular algorithm 
which uniquely determines a particular from a given L, but there does not appear to be any 
reason to prefer this algorithm over other possible ones. Thus, it appears most preferable to 
leave the choice of unspecified apart from the restriction of covariance. In fact, there exist 
two independent sources of ambiguity in 0: 



• As noted above, eq.(|20D allows the freedom to alter by addition of 
an exact (n — l)-form 

+ dY(<f>, 8<j>) (40) 
where the covariant (n — 2)-form, Y, is linear in the varied fields. 

• If we alter the lagrangian by addition of an exact n-form 

L -» L + dfi, (41) 

then the equations of motion are unaffected, so we do not alter the 
dynamical content of the theory. Nevertheless, must be shifted by 

-> + 5fi (42) 

(If is defined by the algorithm implicit in the proof of the above 
lemma, then an additional exact term also would be added to 0.) 

Thus, is ambiguous up to the addition of two terms, 

0^0 + 5/i + dY(4>, 8<f>) (43) 

The consequences of this ambiguity in for the Noether current and charge will be analyzed 
in the next section. 

We conclude this section, by briefly reviewing the definition of the symplectic form, f2, 
in globally hyperbolic spacetimes and investigating its possible ambiguities for asymptotically 
flat solutions. This is of relevance here because f2 is used to define the notion of a Hamilto- 
nian, which, in turn, gives rise to the notions of total energy and angular momentum. Thus, 
ambiguities in f2 could result in ambiguities in these notions. 

Recall that the symplectic current (n — l)-form [14] is defined by 



u(cf>, 5 x <t>, <5 2 </>) = 5 2 0(0, M) - Si&(4>, 5 2 <P) (44) 

Let C be a Cauchy surface. We take the orientation of C to be given by e 01 ... 0n _ 1 = n b ei, ai ^ Mn l 
where n a is the future pointing normal to C and e&ai...a„_i is the positively oriented spacetime 
volume form. We define the symplectic form relative to C by 

fiOM^,^) = / (45) 
Jc 



(More precisely, eq.(45) defines a "pre-symplectic form" on field configuration space. As ex- 
plained in detail in [14], the phase space then is obtained by factoring out by the degeneracy 
submanifolds of f2, and Q then gives rise to a symplectic form on phase space.) If C is noncom- 
pact - as we assume here - then some "asymptotic flatness" conditions must be imposed upon 
the dynamical fields, <f>, (and, hence, on their variations) in order to assure convergence of the 
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integral appearing in (|45|). One normally assumes that the metric, g a b, approaches a flat metric, 
rfab, and the matter fields, ip, approach zero at some suitable rate. The precise asymptotic con- 
ditions appropriate for a given theory depend upon the details of the theory, and, thus, must be 
examined on a case- by-case basis, subject to the following general guidelines: the asymptotic 
fall-off rates of the dynamical fields should be sufficiently rapid that quantities of interest (like 
f2, energy, and angular momemtum) be well defined, but not so rapid that a sufficiently wide 
class of solutions fails to exist. We shall not investigate this issue further here, but will merely 
assume that such suitable conditions have been imposed. 

In principle, the definition of $7 depends upon the choice of C . However, since du = 
whenever 5i<f> and e)2</> satisfy the linearized equations of motion [14], the dependence of on C 



when the equations of motion are imposed is given by an integral of ui over a timelike surface 
near spatial infinity. If sufficiently strong asymptotic conditions at spatial infinity have been 
imposed on the dynamical fields to assure convergence of the integral appearing in (fioT), then 
the integral of u on this timelike surface typically will vanish, and, thus, the definition of Q 
should be independent of C. Of course, if this were not automatically the case, one presumably 
would strengthen the asymptotic conditions imposed upon the dynamical fields in order to make 
Q be independent of C. 

As noted above, is ambiguous up to the terms given in (|43|). The term involving ft does 
not contribute to u or £1, so we find that the only ambiguity in £1 is 

Q -» O + An (46) 

with 

AO= / tfiYOMa^) - <f 2 Y(<Mi0)- (47) 

J oo 

where the integral is taken over an (n — 2)-sphere at spatial infinity. It appears that the 
asymptotic conditions on the dynamical fields needed to ensure the vanishing of Af2 typically 
will be weaker than the conditions needed to ensure that 17, eq.(f45|) is well defined for a given 
choice of 0. In particular, taking account of the difficulties in constructing a covariant (n — 2)- 
form, Y out of the metric and its first variation (as well as e), we see that for a theory in 
spacetime dimension n in which no matter fields are present, the asymptotic conditions 

(n-3) 

9ab ~ Vab + o(r 2 ) (48) 

together with faster fall-off conditions on derivatives of the metric, suffice to ensure that AS7 = 0. 
Thus, it does not appear that the ambiguity in will typically give rise to an ambiguity in the 
definition of 0, for suitable asymptotic conditions on the dynamical fields. 



4 The form of the Noether charge 

In this section, we will obtain an expression for the general structure of the Noether charge 
(n — 2)-form, Q, for a diffeomorphism invariant theory. We begin by reviewing the construction 
of the Noether charge given in || (see also fill ). 

Let £ a be any smooth vector field on the spacetime manifold, M, (i.e., £ a is the infinitesimal 
generator of a diffeomorphism) and let <j> be any field configuration. (</> is not required, at this 
stage, to be a solution of the equations of motion). We associate to £ a and <f> a Noether current 
(n — l)-form, defined by 

J = 0(0,£^)-£-L (49) 
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where 0(0, £^0) denotes the expression obtained by replacing 8(f) with C^4> m the expression 
for 0, and the "centered dot" denotes the contraction of the vector field £ a into the first index 
of the differential form L. A standard calculation (see, e.g., [14|]) gives 

dJ = -E£ ? </> (50) 

which shows J is closed (for all £ a ) when the equations of motion are satisfied. Consequently 
[13] there is a Q locally constructed from cf> and £ a such that whenever <f> satisfies the equations 



of motion, E = 0, we have 

J = dQ (51) 



We refer to Q as the Noether charge (n — 2)-form. Note that for a given J, eq.(pl|) determines 
Q uniquely up to the addition of a closed (and, hence, exact |0|) (n — 2)-form. 

Proposition 4.1 The Noether charge (n — 2) -form can always be expressed in the form 

Q = W c (0)£ c + X cd (0) V [c £ d ] + Y(0, Ctf) + dZ(0, (52) 

where W c ,X. ab ,Y, and Z are covariant quantities which are locally constructed from the indi- 
cated fields and their derivatives (with Y linear in C^(j> an d Z linear in ^). This decomposition 
of Q is not unique in the sense that there are many different ways of writing Q in the form 
$5^), i.e., W C; X ab , Y, and Z are not uniquely determined by Q (see below). However, X afe 
may be chosen to be 

(X c )c 3 ... Cn = —Er c e a bc 3 ...c n (53) 
where E$ cd was defined by eg. (f?ilj, and we may choose Y = Z = 0. 



Proof. We proceed by calculating Q using the choice of given in lemma |3,l| , and the algorithm 
for calculating Q from J given in lemma 1 of [jljj. For the choice of given in lemma 3.1, we 
have 

J = 2E^V d (V 6 6 + V c £ b ) + ®'{<t>, Ctf) — £ • L (54) 



Now the algorithm of lemma 1 of [|13| for obtaining Q from J reduces the highest number of 
derivatives of ^ a appearing in the expression for J by one. Since 0' is linear in the quan- 
tities (Sg a b, 8R a bcd, 'Rabcd, •••) and does not contain any terms involving derivatives of these 
quantities, it follows that ®'(4>,C^4>) is linear in (£ a , V;,£; a ), i.e., it has no dependence on any 
derivatives of £ a higher than first. Since E^ d is antisymmetric in c and d, the term E^V^Vc^f, 
has no dependence on derivatives of £ a . Thus, no derivatives of £ a higher than second appear 



in eq.(5_4), and only the term E^VdV;,^ involves second derivatives of £ a . The contribution of 
this latter term to Q is readily computed, and we find that with our choice of and algorithm 
for calculating Q, we have 

Q = W c (0)e c + X cd V [c ^] (55) 

where W is a covariant (re — 2)-form locally constructed out of the dynamical fields, 4>, and their 
derivatives, and where 

(X cd ) C3 ... Cn = -Eft Cd e abC3 ^. Cn . (56) 

Equation (|55|) gives the general form of Q for our particular algorithm for choosing and 
obtaining Q from J. Recall that had two ambiguities (|43|), one arising from the ambiguity 



in L and the other from its defining equation (]20|) . Using the identity 

£.^ = ^dfi + a?(£ • ft) (57) 
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we see that the ambiguity in gives rise to the following ambiguity in J: 



J-> J + d(f ■/i)+</Y(0,£ c 0) (58) 

Taking into account the additional ambiguity of addition of an exact form to Q, we obtain the 
following ambiguity in Q ||: 

Q Q + £/x + Y(0,£ s 0) + dZ (59) 
Thus, for any choice of Q, we have 

Q = W c f + X C£ V [c ^ + Y(ct>, C^) + dZ. (60) 

as we desired to show. □ 

As stated above, the decomposition ( |52]) of Q is not unique. For example, it is clear that 
for any choice of (n — 2)-form U c (</>), the quantity d(U c (<j))^ c ) can be written as a sum of terms 
of the same form as the first three terms on the right side of (p2;), since we can write it as a 
sum of a term linear in £ c , a term linear in V[ c £</], and a term linear in 2V( C ^) = ££9cd- Thus, 
we can always add the term U c (</>)£ c to Z and make compensating changes in W, X, and Y 
without affecting Q. One might be tempted to impose additional conditions to determine the 
terms W, X, Y, and Z in eq. (|52]) . In particular, it might appear natural to fix the term X ab 
(which plays a key role in the definition of black hole entropy below) by simply requiring it to 
be given by eq.(^) above. However, this proposal suffers from the difficulty that a change of 



Lagrangian of the form (41) - which should have no affect upon the physical content of the 



theory - would, in general produce in a change in X . For this reason, we shall not attempt 



to give unique definitions of the individual terms in eq.(52) but will derive the first law of black 
hole mechanics and give a proposal for defining the entropy of dynamical black holes based only 
upon the general form of Q given in proposition ( |4.1|) . 



5 Examples of Lagrangians and associated Noether currents 
and charges 

In this section, we shall give the symplectic potential 0, the Noether current J, and the Noether 
charge Q arising from three Lagrangians of interest. In giving these examples, we shall simply 
make convenient choices of 0, J, and Q, but, of course, it should be kept in mind that the 
ambiguities (|43|), (pq) , and (|59| ) remain present. 

Our first example is general relativity. We have the Lagrangian 4-form 

L a fecd = TT-^abcdR (61) 
16-7T 

This yields a symplectic potential 3-form 

®abc = e dabc r^g de g fh (V f 5g eh - V e 5g fh ) . (62) 
From this, we obtain the Noether current 3-form 

Jabc = ^e dabc V e (V e £ d] ) , (63) 
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which yields the Noether charge 2-form 

Qab = —^re abcd V c i d . (64) 

1D7T 

Our second example is 2-dimensional dilaton gravity (in the form given in ||), with scalar 
field (ft, coupling constant A, and an additional "tachyon field" T. The Lagrangian 2-form is 

L afe = \e ah e* (R + (Vcft) 2 - (VT) 2 + fi 2 T 2 + a) (65) 
This yields the symplectic potential 1-form 



©a = e ab e 4 



(V b <ft)8<ft - (V b T)5T + ~g bc (V d (6 gcd ) - g de V c (5g de ) - (V d <t>)5g cd + (V c <j>)g de 5g de ) 

(66) 



From this, we obtain the Noether current 1-form 

J a = e ab V c (e^V^ + 2£ [C VV) , (67) 
which yields the Noether charge 0-form (i.e., function) 

Q = - X -e ah (e*V a t b + 2Z a V b e*) (68) 

As our final example, we consider the special case of Lovelock gravity in n dimensions 
obtained by keeping only the terms in the Lagrangian up to quadratic order in the curvature 
(see [17 1). The Lagrangian n-form is 

L ai ...„ n = e ai ... an (j^-R + a(R abcd R abcd - 4R ab R ab + R 2 )^j (69) 

This yields a symplectic current (n — l)-form 

= e dai ,,, an _ 1 {[-^- + 2aR)g de gf h {Vf5g eh -V e 5g fh ) 

+a (-2{V e R)g d ^g ef + 4R de (V e 5g fh )gf h + AR e f(V d 5g ef ) 
-8R ef (V e 6g fh )g dh - 4(V e R d f)5g ef + 4R^ h V h Sg ef )) (70) 



The corresponding Noether current (n — l)-form is 

J ai ...a n ^ = e dai ... an ^V e + AaR)V^C d] + 16a(V f i [e )R d]S + 4aR e <# h V f ^ (71) 
which yields the Noether charge (n — 2)-form 

Q ai ...a n _ 2 = -e deai ... an _ 2 (^V d r + 2a(RV d e + 4V^i? e / + R defh V f ^ (72) 
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6 The first law of black hole mechanics 



In this section, we will use lemma 3.1 and proposition 4.1 to improve upon the derivation of 
the first law of black hole mechanics given in ||. We thereby will prove that the first law 
of black hole mechanics holds for nonstationary perturbations of a black hole in an arbitrary 
diffeomorphism covariant theory of gravity, without any restriction on the number of derivatives 
of fields which appear in the Lagrangian. 

Let be any solution of the equations of motion, and let 50 be any variation of the dynamical 
fields (not necessarily satisfying the linearized equations of motion) about 0. Let £ a be an 
arbitrary, fixed vector field on M. We then have || 

53 = 50(0,£ 5 0)-£-5L 

= 50(0,££0) -£-d0(0,50) 

= 50(0,££0)-£ c 0(0,50) + d(£-0(0,50)) (73) 



where eq.(|2Cj) together with E = was used in the second line, and the identity (|57| ) on Lie 
derivatives of forms was used in the last line. Since our choice of is covariant, £^0 is the 
same as the variation induced in by the field variation 5'0 = £^0. Consequently, we have 

50(0, £ c 0) - £ 5 0(0, 5<j>) = w(0, 50, £ c 0) (74) 

where u was defined by eq. (fH|) . We therefore obtain 

u>(0,50,£ 5 0) = 5J-d(£-0) (75) 

The fundamental identity which gives rise to the first law of black hole mechanics applies to 
the case where £ a is a symmetry of all of the dynamical fields - i.e, £^0 = - and 50 satisfies 
the linearized equations of motion. When £^0 = 0, the left side of eq.(|75|) vanishes, and when 
50 satisfies the linearized equations, we may replace 5J by 5dQ = d5Q on the right side. Thus, 
we obtain, 

d5Q - d(f • 0) = (76) 
Integrating this equation over a hypersurface, S, we obtain 

/ $Q[£]-£. 0(0,50) =0 (77) 

Jd~ 

We emphasize that the only conditions needed for the validity of eq. ( |77|) are that be a solution 
to the equations of motion, E = 0, satisfying £^0 = 0, and 50 be a solution of the linearized 
equations (not necessarily satisfying £^50 = 0). 

We shall be interested here in the case where H is an asymptotically flat hypersurface in an 
asymptotically flat spacetime. In this case, a boundary term from an asymptotic (n — 2)-sphere 
at infinity will contribute to eq.([77[). The following argument shows that this boundary term 
has the natural interpretation of being the variation of the "conserved quantity" canonically 
conjugate to the asymptotic symmetry generated by £ a . 

Consider a solution, 0, corresponding to an asymptotically flat, globally hyperbolic space- 
time, with Cauchy surface, C, having a single asymptotic region and a compact interior. We 
return to eq.([75|) but no longer impose the additional assumptions that £^0 = or that 50 sat- 



isfy the linearized equations of motion. We integrate eq.(75) over C taking into account eq.fliq) 



14 



and the fact that, by definition, Hamilton's equations of motion for the dynamics generated by 
the time evolution vector field £ a are 

5H = n((j),6(t),Cz<f>) (78) 
We thereby find that if a Hamiltonian, H, exists for the dynamics generated by £ a , then 

5H = 6 [ J - [ d(£-&) 
Jc Jc 



S [ J - / £ • (79) 

JC Joo 



Thus, a Hamiltonian for the dynamics generated by £ a does exist if (and only if) we can find a 
(not necessarily diffeomorphism covariant) (n — l)-form, B, such that 



5 I £ • B = / £ (80) 

J oo J oo 



in which case H is given by 



H= [ J - / £ B (81) 

JC Joe 

Now evaluate H on solutions. We then may replace J by (iQ, whence H becomes 

H = [ (Q - £ • B) (82) 



Thus, we have shown that in any theory arising from a diffeomorphism covariant Lagrangian, 
the Hamiltonian - if it exists - always is a pure "surface term" when evaluated "on shell". 
Similarly, for a closed universe (i.e., compact C), the Hamiltonian always vanishes "on shell". 

We now shall assume that the asymptotic conditions on the dynamical fields have been 
specified in such a way that when £ a is an asymptotic time translation, B exists, and the 
surface integrals appearing in eq.(^2|) approach a finite limit at infinity. We define the canonical 
energy, £ to be the value of the Hamiltonian, i.e., 

£ = [ (Q[t]-t-B) (83) 



where t a is an asymptotic time translation. We then adopt eq.(|83|) as the definition of the 
canonical energy associated to any asymptotically flat region of any solution, whether or not 
the spacetime is globally hyperbolic. 

We illustrate this definition of canonical energy by evaluating £ for vacuum general relativity. 
We consider spacetimes which are asymptotically flat in the sense that there exists a flat metric 
rjaj, such that in a global inertial coordinate system of n a b we have 



g^u = v + °0-/ r ) (84) 

and 

0(l/r 2 ). (85) 



dx c 

Let t a be the asymptotic time translation (d/dt) a , and let the 2-sphere at infinity be the limit as 
r — > oo of the coordinate spheres r,t = const. Then from our previously calculated expression 
for Q a b, eq.(64), we find that 

/ Q[i] = / e abcd VH d 

Joo 107T Joq 

1 W 



167T Joo V dr dt 
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Note that for a stationary spacetime with stationary Killing field t a , the first line of eq. ( pq ) 
shows that Q[t] is precisely one-half of the Komar mass (see, e.g., fL2| ). 

We now compute the contribution to £ from the second term on the right side of eq.(83). 



Using eq . (|62|) , we have 

Jj a ®abc = J^dSr d g de gf h (X7 f 5g eh -V e 5g fh ) 

= J dSg rr (g tt {d t 5g rt - d r 8g tt ) + h l ^(d t 6h r] - d r bh tJ \ 

= -^6 J dS ((d r9tt - d t g rt ) + r k 0{dih kj - dkhj)) (87) 

where r a = (d/dr) a and fay is the spatial metric. Thus, we see that eq.(|S0|) holds if B is chosen 
to be any 3-form such that asymptotically at infinity, we have 

t a B abc = -r^hc ((dr9tt ~ dtgrt) + r k h ij (dih kj - dkhij)) (88) 

where ef, c is the volume 2-form for the sphere at infinity. Combining this with (^) we find the 
that the canonical energy, £, for general relativity is 

£ = [ Q[t]-t-B 

J oo 

= T7T / dSr k 0(dih kj - 8^) 

l07T Jog 

= M ADM (89) 

where denotes the ADM mass. Thus, the term t ■ B cancels the contribution to £ from 

the term Q, and, in addition, provides the term Madm-, thereby making our definition of £ in 
vacuum general relativity reduce to the standard, ADM, definition of energy. Note, however, 
that additional contributions to £ in general relativity can occur when long range matter fields 
are present; see [jlj] for an explicit evaluation of the contribution to £ for Yang- Mills fields. 

When £ a is an asymptotic rotation, ip a , we may choose the surface at infinity to be every- 
where tangent to ip a , in which case the pullback of ip ■ © to that surface vanishes. Hence, we 
define the canonical angular momentum, J of any asymptotic region by 

J = - I QM (90) 



where it is assumed that the asymptotic conditions on the dynamical fields are such that this 
surface integral approaches a well defined limit at infinity. (The relative sign difference occuring 
in the definitions ( |8^ ) and (|90| ) traces its origin to the Lorentz signature of the spacetime 
metric. The same relative sign difference occurs in the definitions, E = —p a t a and J = +p a ip a , 
of the energy and angular momentum of a particle in special relativity.) In the axisymmetric 
case in vacuum general relativity, eq. ( |90| ) is precisely the Komar formula for angular momentum. 
Thus, we see that in any theory, the Komar- type expression — Q[(p] always yields the angular 
momentum, but / Q[i] does not, in general, yield the energy. Indeed, since the Komar and 



ADM masses agree for stationary solutions in general relativity [18], we see from our calculation 
above that f Q[t] yields only half of the energy in that case. It is the presence of the "extra" 
t ■ B term in eq. ( |83|) which accounts for this well known "factor of 2" discrepancy in the Komar 
formulas for mass and angular momentum in general relativity. 
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We now are ready to apply eq.(77) to the case of a stationary black hole solution with 
bifurcate Killing horizon. Let £ a be the killing field which vanishes on the bifurcation (n — 2)- 
surface E, normalized so that 

e=t a +^^) (9i) 

where t a is the stationary Killing field with unit norm at infinity, and summation over [i is 
understood. (This equation picks out a family of axial Killing fields, acting in orthogonal 

planes, and also defines the "angular velocities of the horizon", . No summation is required 
when the spacetime dimension is less than five, and, of course, the the second term on the right 
side is entirely absent in two dimensions.) Let H be an asymptotically flat hypersurface having 
E as its only "interior boundary". Then, taking into account eq.(|9~l|), the definitions of £ and 
J , and the fact that £ vanishes on E, we obtain directly from eq.([77|) the result 

5^Q[C]=5£-nP5J M . (92) 

We now are ready to state and prove the first law of black hole mechanics in a form which 
strengthens the results of || by establishing the general validity of this law for nonstationary 
perturbations. 

Theorem 6.1 Let (j) be an asymptotically flat stationary black hole solution with a bifurcate 
killing horizon, and let 5(f) be a (not necessarily stationary), asymptotically flat solution of the 
linearized equations about <f>. Define S by, 

S = 2tt j X cd e cd (93) 



where ~K cd is as given in proposition 4-1, and the integral is taken over the bifurcation (n — 2)- 
surface, E, with e C( i denoting the binormal to E (i.e., e is the natural volume element on the 
tangent space perpendicular to E, oriented so that e c dT c R d > when T a is a future- directed 
timelike vector and the spacelike vector R a points "towards infinity"). Then we have 



K 



5S = 5£- n ( ^5J M (94) 



2tt 

where k is the surface gravity of the black hole. 



Proof. The theorem will follow from eq.(|92|) provided that we can show that 

6 L m= ^ 6s (95) 



To evaluate the left side of this equation, we appeal to proposition |4.1| and examine the contri- 
bution of each of the four terms individually. Since £ vanishes on S, it is clear that the term 
W c £ c contributes neither to Q nor to its variation. Similarly, the term dZ clearly also makes 



no contribution to the left side of eq.(95). Since C^(f) = 0, the term Y vanishes in the stationary 



background, and its first variation is given by 



SYfaCtf) = Y{^C t 5cj>) 
= CpYU 



£ ■ dY + d(£ • Y) (96) 
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where the Lie derivative identity (|57]) was used in the last line. It follows immediately that the 
term Y also makes no contribution to the left side of eg. (|9~5|) . Thus, we have 



5 / S Q[£]=^/ 2 XC ^) V [^] (97) 

Now, in the stationary background, we have, on S 

V c £ d = Ke cd (98) 

Furthermore, since £ a = on E, and <5£ a = everywhere, we have 

5V c £ d = (99) 

on E. Consider, now, the variation, 5e c d of the binormal, e c d , with an index raised. Clearly, 
s c 5e c d = for all s c tangent to E, so 5e c d has no "tangential-tangential" piece. However, since 
e c d ed c does not vary as the metric is changed, it follows that g a [c^ e d] a nas no "normal-normal" 
piece with respect to the background metric. Thus, writing 

w cd = ^[c^d] ~ Ke cd (100) 
we have that w cd vanishes in the stationary background, and 

Swd = $ (g a [d(V c] f - Ke c] a )) 

= -Kg a [ d Se c] a (101) 

so that Sw cd has only a "normal-tangential" piece with respect to the background metric. Thus, 
substituting in eq. (|9~7|) , we find, 

6 [ Q[f] = s [ ^ cd m^ cd +w cd ] 

= ^5S+ [ X cd 5w cd (102) 



Finally, we note that since C^cj) = 0, we have C^X. cd = 0, and, hence, by lemma 2.3 of [19|, at 
each point of S, X cc( must be invariant under "reflections" about S, i.e., X crf must be invariant 
under the map of the tangent space which reverses the normal directions to S but keeps the 
tangential directions unchanged. On the other hand, since 5w cd is purely "normal-tangential", 
it reverses sign under reflections about E. However, the pull-back of X. cd 5w cd to E is purely 
tangential, and, hence, invariant under reflections. Consequently, the pull-back of X. cd 5w cd to 
E must vanish, so the second term on the right side of eq. (|102| ) does not contribute. □ 

It should be noted that in the above discussion, E was explicitly chosen to be the bifurcation 
surface of a bifurcate Killing horizon. However, as pointed out in for a stationary black hole 
with bifurcate horizon, the integral of Q is independent of the choice of horizon cross-section. 
Namely, the difference between the integrals of Q over cross-sections E and E' is given by an 



integral of J over the intervening portion of the horizon. However, by eq.(4£), the pullback of 
J to the horizon vanishes, since C^cf) = and the pullback of £ ■ L vanishes since £ a is tangent 
to the horizon. 

Furthermore, if we define the entropy, S, for an arbitrary horizon cross-section, E', of a 
stationary black hole by 

S[E'] = 2tt/ X cd e' cd (103) 
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where e' cd denotes the binormal to then S also is independent of the choice of S' ||. To 
prove this, we note that since X cd is invariant under the one parameter group of isometries, Xt, 
generated by £ a , it follows immediately that S[xt{^')] = However, as t — > — oo, Xt(^0 

continuously approaches the bifurcation surface, S, and (since ~K cd is smooth) we thus obtain 
S[T,'] = S[Ti], as we desired to show. It follows immediately that for stationary perturbations, 
the first law of black hole mechanics fl94| ) holds with S taken to be the entropy of an arbitrary 
cross-section of the horizon. However, when nonstationary perturbations are considered, it is 



essential for the validity of eq.(94) that S be evaluated on the bifurcation surface, S. 

As emphasized at the end of the section 4, the decomposition of Q given by eq.([52"[) does 
not uniquely determine X crf . Nevertheless, theorem 6.1 and its proof show that all of the 



different possible choices of X crf yield the same value of the entropy, S, for a stationary black 
hole. Furthermore, even for nonstationary perturbations, the first variation, 5S, of S on £ is 
independent of the choice of X cc( . However, for nonstationary perturbations, SS will, in general, 
depend upon the choice of X crf when evaluated on an arbitrary cross-section, of the horizon, 
and the dependence of S upon the choice of X crf becomes even more severe if we attempt to 
generalize the notion of entropy to an arbitrary cross-section of a nonstationary black hole via 
eq.(|103|). We turn, now to an analysis of the definition of entropy for nonstationary black holes. 



7 A prescription for dynamical black hole entropy 

In this section we will suggest a definition of the entropy, S dyn , for a "dynamical" (i.e., nonsta- 
tionary) black hole. We seek a formula of the general type 

S dyn [C] = I ± cd ^)e cd (104) 



c 

where C is an arbitrary cross-section of the event horizon of a dynamical black hole, and X crf is 
a diffeomorphism covariant (n — 2)-form locally constructed out of the dynamical fields, 4>, and 
their derivatives by an algorithm whose sole input is the Lagrangian, L. There are four basic 
criteria which our definition of Sd yn must satisfy: 

1. For an arbitrary cross-section, £', of a stationary black hole, we must 
have 

S dyn [Z'} =5[E'] =2tt f X cd e' cd (105) 
Jt,' 



(see eq.(103) above). 



2. For an arbitrary (nonstationary) perturbation of a stationary black 
hole, on the bifurcation surface, E, we must have 

5S dyn m = 5S = 2<KdJ^X cd e cd (106) 

(see eq.(|93|) above). 

3. If we alter the Lagrangian by addition of an exact n-form 

L -> L + dfx (107) 

then the definition of S dyn should not change, since there is no change 
in the dynamical content of the theory. 
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4. At least for an appropriate class of theories, Sdyn should obey a "second 
law", i.e., Sdyn should be a non-decreasing quantity when evaluated on 
successively "later" cross-sections of the horizon of a dynamical black 
hole. 

The last of these criteria is by far the most interesting and important. Unfortunately, it 
also is the most difficult to analyze in a general theory of gravity for at least the following two 
reasons: First, it seems clear that, unlike the "first law", any proof of the second law would 
need to make detailed use of the equations of motion of the theory. Second, it seems clear 
that the "second law" should hold only for the case of theories which satisfy certain physically 
reasonable criteria, likely examples of which are the existence of a well posed initial value 
formulation, cosmic censorship, and the property of having positive total energy. For example, 
even for general relativity, the second law can fail if matter is present which fails to satisfy the 
weak energy condition. However, it is far from clear as to precisely what conditions should be 
imposed upon a theory for the validity of the second law to hold, and - even if these conditions 
were known - it undoubtedly would be highly nontrivial to determine whether a given theory 
satisfied them. 

Despite these two difficulties, there are some hints that it may be possible to prove some 
general results pertaining to the second law. In particular, we saw in the previous section that 
the entropy, S, of a stationary black hole is just its Noether charge with respect to the horizon 
Killing field, £ a . Thus, the change in entropy between cross-sections C and C of a stationary 
black hole is given by the flux of the corresponding Noether current through the horizon between 
C and C For a stationary black hole, this flux, of course, vanishes. However, if Sdyn could 
similarly be identified as the Noether charge of an appropriate vector field, one might be able to 
establish a relationship between the "second law" and positive energy (i.e., positive net Noether 
flux) properties of the theory. Another suggestive fact is that the quantity X cc( which plays a 
key role in the definition of entropy for stationary black holes can be chosen to be very simply 
related to E R bcd (see eq.([53|) above), and E^ > , in turn, is a term in the equations of motion 
(see eq,(|37|) above). Thus, there is a hint that it may be possible to define Sd yn in such a way 
that its dynamical properties may be directly related to the equations of motion of the theory. 
Unfortunately, we have not, as yet, succeeded in developing either of these hints into any results 
regarding proposed definitions of Sdyn- Thus, for the remainder of this section, we shall not 
consider criterion (4) further, and will merely seek a definition of Sd yn which satisfies conditions 
(l)-(3). 

An obvious first try at defining Sd yn via an equation of the form ( |104D would be to simply set 
X cd = X crf , with X cd given by the decomposition (^) of Q. However, we already emphasized 
above that this decomposition is not unique. Although - as discussed at the end of the previous 
section - this ambiguity does not affect the evaluation of S on an arbitary cross-section of a 
stationary black hole or the evaluation of 5S on the bifurcation surface of a stationary black 
hole, this ambiguity in ~K cd is of importance for a dynamical black hole. 

An obvious try at circumventing this difficulty would be to continue to set ~K cd = X crf and 
simply fix ~K cd by some definite algorithm. In particular, the choice 

Xcd Tpabcd r (A(\q\ 

a 3 ...a n — e aba 3 ...a n - 

(see eq.([53"l) above) appears to be particularly simple and natural. This proposed definition of 
Sdyn clearly satisfies conditions (1) and (2) above. However, it is not difficult to verify that it 
fails 0] to satisfy condition (3): By adding an exact form to L which has suitable dependence 
upon the curvature, we can alter E R bcd in such a way as to produce nonvanishing changes in 
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Sdyn for nonstationary black holes. We feel that it is unlikely that any other simple algorithm 
for fixing ~K cd for a given L will fare any better in this regard. 

Thus, it is a nontrivial challenge to find any prescription for Sd yn of the form Q104| ) which 
satisfies conditions (l)-(3). We now shall demonstrate that such a prescription does exist. The 
basic idea will be to construct new dynamical fields relative to a cross section C, which make C 
"look like" a bifurcation surface of a stationary black hole. We then shall define Sd yn [C\ to be 
the entropy of this stationary black hole. Before giving a precise statement of our prescription, 
we give the following two definitions: 

Definition 7.1 ; Let C be a (n — 2) -dimensional spacelike surface in an n- dimensional space- 
time, and let M ai "' ak b 1 ...bi be a (spacetime) tensor field defined on C. Then M ai "' ak i )1 ,..i >l is said 
to be boost invariant on C if, for each p G C, M ai ~ ak b 1 ...b l is invariant under Lorentz boosts in 
the tangent space at p in the 2- dimensional timelike plane orthogonal to C. 

The following simple criterion can be used to check if a tensor field M ai " afc ; )1 ^ ( is boost 
invariant on C. At each point p £ C, choose a null tetrad with null vectors l a and n a orthogonal 
to C, and spacelike vectors s® tangent to C. Expand M ai, ^ ak b 1 ...b l in this basis. Then it is easy 
to verify that iW ai " , ° fc & 1 ,..b, is boost invariant if and only if its basis expansion is "balanced" 
with respect to l a and n a , i.e., if the basis expansion coefficients are nonvanishing only for terms 
involving equal numbers of l a, s and n a 's. This motivates the following definition. 

Definition 7.2 ; Let C be a (n — 2) -dimensional spacelike surface in an n-dimensional space- 
time, and let M ai, " ak b 1 ...b l be a (spacetime) tensor field defined on C. We define the boost 
invariant part of M ai " Mk b 1 ...bi to be the tensor field on C obtained by keeping only the terms 
which are balanced with respect to l a and n a in a null tetrad basis expansion. 

It is easily seen that the boost invariant part of M ai '" ak b 1 ...b, does not depend upon the 
choice of null tetrad appearing in the definition. 

Note that the spacetime metric, g a b, on C is automatically boost invariant. However, the 
curvature of g a b and its derivatives need not be. Nevertheless, we may define a notion of the 
boost invariant part (up to order q), gl, of the spacetime metric in a neighborhood of C. The 
curvature of g^ b and its covariant derivatives up to order (q — 2) then will automatically be 
boost invariant on C. This construction of g | will lead directly to a proposal for defining Sd yn . 

To define g 1 ^, it is convenient to introduce a coordinate system in a neighborhood of C as 
follows [19[| . Define a null tetrad l a ,n a ,s^ on C as above, with l a n a = —1. Let O be any 
neighborhood of C sufficiently small that each point ieO lies on a unique geodesic orthogonal 
to C. Given x € O we find the point p € C and the geodesic tangent v a in the (n — 2)-plane 
normal to C such that x lies at unit affine parameter along the geodesic determined by p and 
v a . We assign the coordinates (U, V, s±, ...s n -2) to x G O by taking (U, V) to be the components 
of v a along l a and n a , respectively, and taking Sj to be (arbitrarily chosen) coordinates of p 
on C. We denote by d a the flat derivative operator associated with these coordinates. Note 
that a change in tetrad, l a — > al a , n a — > a~ l n a , at p € C (corresponding to a Lorentz boost 
in the tangent space in the plane orthogonal to C) induces the linear change in coordinates, 
U — ► a~ l U , V — > aV, Si — ► Si. Since linearly related coordinate systems define the same 
"ordinary derivative operator", it follows that d a does not depend upon the choice of l a and n a , 
and so is invariant under the action of Lorentz boosts in the plane orthogonal to C. 
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Now consider the first q terms in the Taylor series expansion of g ab around C in U and V 



a® (a?)- T ^ V &m+ngaP (s ) 

y<lb I mini ^ Z>rnTTXn-lA ai ) 



n,m=0 



ail 



d m Ud n V 



(dx a ) a (dxP) b . 



(109) 



u=v=o 



The coefficients appearing in this expansion are just components of the tensors d Cl ...d Cr g a b on 
C, namely, 



/ I am. T r an. T r \ 1 ' 



a/3 



gmjjQny 



(dx a ) a (dx?) b = l c 



u ci-- u c m+n 



9ab 
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u=v=o 



We define g^ b by replacing each tensor, d Cl -.-d Cr g a b, appearing in the expansion of g^ b by its 
boost invariant part. In other words, we alter g ab by extracting the boost invariant part of the 
coefficients of the first q terms of its Taylor expansion in U and V. 

The nature of g^ b can be best elucidated in the case where g a b is analytic, in which case we 
may set q = oo and write g^ b for g 1 ^ . It then follows that the vector field 



dU' 



dV 



(111) 



(which induces Lorentz boosts of the coordinates) is a Killing field of the metric g T ab , that is 



o. 



(112) 



Furthermore, ^ a vanishes on C. Thus, our construction of g^ b has, in effect, created a new 
spacetime (which is not necessarily a solution of the field equations) in which C is the bifurcation 
surface of a bifurcate Killing horizon. 

In an exactly similar manner, we define the boost invariant part, ip Iq of the matter fields tp 
(up to order q) by extracting the boost invariant part of the coefficients of the first q terms of 
the Taylor expansion of t\> in U and V about C. It then follows that £ also Lie derives ip 9 up 
to order q. 

Our proposal for defining Sd yn is the following: Choose q to be larger than the highest 



derivative of any dynamical field appearing in the decomposition of Q given in proposition 4.1. 
Given a cross-section, C, of the horizon of a black hole, we replace g ab by g ab and ip by ip Iq in 
a neighborhood of C. Define Q[£] on C to be the Noether charge (n — 2)-form of the dynamical 
fields cf) Iq = (g a l,4> Iq ) for the vector field £ a defined by eq.( |lll| ) above. Define Sd yn at "time" C 

(113) 



S dyn [C] =2vr / Q[£] 



Equivalently, by proposition 4.1 we have 



Sdyn[C] = 2-K 



X cd e cd 



(114) 



where 

x cd (0) = y, cd {4> lq ) (us) 

Equation ( |114[ ) shows that Sdyn 15 °f the desired general form (104), and the equivalence of 
eqs. (|ll3p and flH4|) shows that the right side of (114) does not depend upon the choice of 
X crf in the decomposition of proposition 4.1. Note, incidentally, that since X crf is a nonlinear 
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function of the dynamical fields 4>, the tensor field X crf is not necessarily equal to the "boost 
invariant part" of the tensor field ~K cd (</>). (Use of the boost invariant part of X afe (0) would not 
yield a satisfactory prescription for Sd yn since it would, in general, fail to satisfy condition (3) 
above.) More generally, for a nonlinear tensor function (3 of the dynamical fields eft, we have, in 
general, [(3(cp)] Iq ^ (3{4> Iq )- On the other hand, if [3 is linear in (ft, then [f3(4>)] Iq = f3(4> Iq ). 

We now may verify that our definition of Sd yn satisfies conditions (l)-(3) above. First, if C 
is taken to be the bifurcation surface, S, of a stationary black hole, then (j) Iq = 0, so, clearly, 
(S^nP] = S[Ti]. On the other hand, if £' is an arbitrary cross-section of a stationary black 
hole, then since our prescription for defining Sdyn is a "local, geometrical" one, by isometry 
invariance we clearly have Sd yn [Xti^')\ = Sdy n [^')- But it also is clear that our prescription 
for defining Sd yn is such that Sd yn [C] varies continuously with C. From these facts, it follows 
immediately by the same argument as given below eq.( |103|) that Sd yn [^'] = <5%nP]- Thus, we 
have 

SdyrJtf] = SdynlZ} = S[Z] = 5[S'] (116) 

i.e., condition (1) is satisfied. 

To verify that condition (2) holds, we note that since we have 4> Iq = 4> on the bifurcation 
surface, S, of a stationary black hole, and since <5[X cd e c d] clearly is linear in 5<fi, it follows 
that <5[(X crf — X. cd )e c d] has no boost invariant part. However, this immediately implies that the 
pullback of this differential form to X vanishes, from which it follows that 5Sd yn P] = <^>p]> as 
desired. 

Finally, the complete ambiguity in Q (including that arising from the change in Lagrangian 
L — > L + dfj.) is given by eq. (|59|) . It is manifest that none of these ambiguous terms can 



contribute to J c Q[£]. Consequently, we see from eq. (|ll3[) that condition (3) holds. 

Thus, we have proven the existence of a definition of Sd yn which satisfies conditions (l)-(3). 
These conditions do not uniquely determine Sdyn- Nevertheless, we have been unable to come 
up with any "natural" alternative definitions of Sd yn - Thus, we believe that our definition of 
Sdyn is a serious candidate for the definition of the entropy of a nonstationary black hole in a 
general theory of gravity. 

We conclude this section by evaluating Sdyn for the three theories considered in section 5. 
Consider, first, vacuum general relativity. Let C be an arbitrary cross-section of a black hole, 
let e a b be the binormal to C, and let e a b denote the volume element on C. Comparing eqs.(|5z 
and (64), we see that the 2-form ~K cd is given by simply 



(X cd ) ab = --^ r e ab cd (117) 

Since ~K cd does not depend upon any derivatives of g a b, it is clear that it is unaffected when g a b 
is replaced by its boost invariant part. Thus, we obtain, 



Sdyn [C] — g J e ab e cd 



7 / e - 

4 Jc 

Area[C] 



(118) 



in agreement with the usual formula for the entropy of a dynamical black hole in general 
relativity. By the area theorem, this definition of Sd yn satisfies the "second law" (assuming that 
the cosmic censor hypothesis is valid). Note that if we add to the Lagrangian "matter terms" 
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which have no explicit dependence upon the curvature, then X crf does not change (see eq. 
so eq. (|118| ) also holds for general relativity with matter present, provided only that the matter 
does not have an explicit coupling to the curvature in the Lagrangian. 

The calculation of Sdyn for dilaton gravity with Lagrangian (|65|) in two spacetime dimensions 
proceeds similarly. We see from eq. (|6^) that the 0-form X cd is given by 

X cd = -\^>e cd (119) 

Again ~K cd does not depend upon any derivatives of the dynamical fields, and is unchanged 
when they are replaced by their boost invariant parts. In this case, a cross-section, C, of the 
horizon is a point, and we obtain 

S dyn [C] = 2vr^| c (120) 

It is known that this definition of Sd yn also satisfies the second law j|. 

Lovelock gravity provides a more interesting illustration of our prescription, since it can be 
seen from eq.(72) that X cd contains terms involving the curvature, 



(X_ cd ) — -e cd ( — !— + 2aR I - Hne^ d t R c \ f -2aeth R cdfh (121) 

{sl )ai...a„-2 — fc a\...a n —2 I j^g^j- I out fa 1 ...a n -2 rL ^« fc /ftai...a n _2- rl K 1 ^ 1 ) 

and the replacement of the metric by its boost invariant part will have a nontrivial effect. 
Indeed, since, after this replacement is made, both extrinsic curvatures of C embedded in M 
vanish, we see (using a "Gauss-Codazzi" equation - see, e.g., fl~2]] ) that the curvature of the 
boost invariant part of the metric satisfies, 

= R- 2t ab R ab + t ac t bd R abcd (122) 

where t ab = —n a ni, + r a r b is the metric for the subspace orthogonal to C (spanned by the unit 
timelike and spacelike normals, n a and r a respectively) and ( n ~ 2 ^R is the scalar curvature of C. 



From eqs.( |121| ) and (122), we obtain 



c V.a„ = + 4« ( n - 2 W*]) e ai ...a n _ 2 (123) 



t-cd^*- ui...u n l g 

where ( n-2 ) R[g Iq ] is the (n — 2)-scalar curvature of C computed with the boost invariant part of 
the metric and e ai ... an _ 2 is the volume form C. However, we clearly have ( n ~ 2 ^ R[g Iq ] =( rt_2 ) R[g\. 
Hence, we obtain, 

Sdyn = ^Area[C] + 8ira (n " 2) R (124) 
Note that this formula differs from what would be obtained from simply substituting the ex- 



pression (121) into eq.(^). It is not known whether this definition of Sd yn satisfies the second 
law. 
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Appendix: Applications to Theories with a Nondynamical 
Metric 



In the body of this paper, we have considered theories which are diffeomorphism covariant in the 
sense of eq. (||) . It was seen in section 2 that this condition implies the absence of "nondynamical 
fields" in the Lagrangian. In particular, the diffeomorphism covariance condition excludes 
the case of theories with a nondynamical metric, such as theories of fields in flat spacetime. 
Nevertheless, a number of formulas and results derived in the body of this paper continue to 
hold for theories with a Lagrangian locally constructed out of a metric, g ab , and matter fields 
ip, of the form @, i.e., for the Lagrangian 

L = L (gab,VaiRbcde, ■■; V( 0l ...V^Rbcde, t/>, V ai ^, V( 0l -V ai )tpj (125) 

but where the metric, g ab , is now treated as a fixed, nondynamical entity, so that, in particular, 
the equations of motion, E 9 = 0, no longer are imposed. The purpose of this appendix is 
to present simple, unified derivations of some formulas and results (most of which are "well 
known") for such theories with a nondynamical metric. 

In a theory with Lagrangian of the form ( |125D but with nondynamical metric, we define the 
stress-energy tensor, T ab = T^ ab \ of the matter fields by, 

T ab e = 2(E g ) ab (126) 

For each vector field, we again define the Noether current, J, by eq,(|49|) above. However, 
the (matter) equations of motion no longer imply that J is closed. Indeed, by eq.(^), we see 
that when Ew, = 0, we have 

dJ = -(B g ) ab C^g ab 
= -T ab V {a £ b) e 
= -Va[A]e + aV s [T a6 ]e 

= -d(k-e)+t b V a [T ab ]e (127) 

where 

By inspection of eq.( |127j ), we see that the n-form a [T ab ]e is exact for all £ a . However, since 
£ a is arbitrary, this is impossible unless 

V a T ab = (129) 

which shows that the stress-energy tensor is covariantly conserved whenever the matter equa- 
tions of motion hold. Note that eq. (ll27|) then yields simply 

d(J + k ■ e) = (130) 

from which it follows immediately |]l3| that 

J + k ■ e = dK (131) 

where K is locally constructed out of g ab , t/j, £ a , and their derivatives. In other words, we have 
shown that, apart from a "surface term", the Noether current is equivalent to the stress-energy 
current —T ab ^ b . 
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It is important to note that the equations of motion for g a b were not used anywhere in the 
derivation of eqs. (ff3j) - (l75|) or eqs.(|78|)-(^T|). Thus, these equations remain valid in the case of a 
theory with a nondynamical metric. In particular, if a Hamiltonian exists for a time translation 
vector field, t a , on a globally hyperbolic, asymptotically flat spacetime, then it is natural to 
define the canonical energy at "time" C by 

£ = [ J - / t B (132) 

JC Joe 

(see eqs,(|8l|) and ( |83"| ) above). In other words, apart from the possible "surface term" t ■ B 
(which vanishes in most of the commonly considered theories of matter fields in a background 
spacetime), the canonical energy is simply the integral of the Noether current, J, over a Cauchy 
surface. However, since we no longer have J = dQ, this volume integral no longer can be 
converted into a surface integral. In particular, £ depends upon the choice of t a in the interior 
of the spacetime, not just upon its asymptotic value at infinity. Note also that since J need not 
be closed (see eq.( |127D ), £ need not be conserved, i.e., independent of C. 
Using eq.( 13ip , we find 



/ k-e+ [ (K — i • B) 

JC Joe 



'C 

T ab n a t b e+ f (K-t-B) (133) 

C Joo 

where n a denotes the future-directed unit normal to C, and €b l ...b n _ 1 = n a e a b 1 ...b n _ 1 is the 
natural volume element on C. Thus, apart from some possible surface term contributions which 
can arise from both K and B, the canonical energy is given by the usual formula involving an 
integral of the stress-energy tensor over C. 

As noted above, in general £ is not conserved, i.e., independent of choice of Cauchy surface, 
C . However, if the spacetime metric is stationary, Ctdab = 0, (but stationarity need not be 
imposed upon the matter fields), then the first line of eq, (|127l ) shows that J is closed. Equation 



(131) then immediately implies that the stress-energy current form —k ■ e also is closed (as also 
could easily be verified directly). Equation ( |l32j ) then implies that £ does not change when 
the Cauchy surface, C, undergoes variations of compact support. In the usual case where t ■ B 
vanishes at infinity and J goes to zero suitably rapidly at infinity, £ will take the same value 
for all asymptotically flat Cauchy surfaces. 

Now, suppose that g a b is stationary - i.e., Ctg a b = - and suppose that ip(X) is a one- 
parameter family of solutions to the matter equations of motion (in the fixed metric g a b) such 
that tp(0) is stationary, i.e., £tip(0) = 0. Let £(A) denote the canonical energy of these solutions. 
Then, since Ctdab = 0, by eqs. @) and (H), we have for all A 

^(A) = fiWA),^^(A)) (134) 

In particular, since Cttp vanishes at A = 0, we see that the first variation of the canonical energy 
about a stationary solution vanishes, 

5£ = (135) 

Now take the derivative of eq.( |134D with respect to A and evaluate the resulting equation at 
A = 0. A nonzero contribution will occur on the right side only when the A- derivative acts on 
Ctip. We thereby find that the second variation of canonical energy about a stationary solution 
ip is given by 

1 f] 2 £ 1 

s 2 £ = 2^U=o = -nftMV,A&/>) (136) 
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where 

W = £l« (137) 

Note, in particular, that 5 2 £ depends only upon Sip, and not upon 8 2 ip. 

Equation (|136| ) is one of the key results of this Appendix. To elucidate its meaning, we 
note that if Sip satisfies the linearized equations of motion about a stationary solution, ip, then 
so does C t Sip. Consequently, the symplectic current form u(ip,5ip, C t Sip), defined above by 



eq.(44), is closed. Thus, its integral over a Cauchy surface, C, yields a conserved quantity for 



perturbations. Equation (136) shows that, apart from a factor of 2, this conserved quantity 
is just the second order change in the canonical energy associated with this perturbation. By 
our previous results, we see that this conserved quantity is equivalent - up to possible "surface 
terms" - to the conserved quantities J c S 2 J and j c S 2 T a t>n a t b e. 

As a simple application of the above result, consider the theory of a linear field ip in a 
stationary spacetime, where by "linear" we mean that L is quadratic in ip, so that the equations 
of motion for ip are linear. In this case, the equations of motion are the same as the linearized 
equations about ip = 0, so we may choose the "unperturbed solution" to be ip = 0, and we may 
write Sip = ip in the above formulas. We also have 5 2 £ = £ and 5 2 T a b = T^. Hence, we obtain 
from eqs. (|l33|) and (|136[) 

Q{Sip,C t 5ip) = 2 ( T ab n a t b l + 2 f (K - t ■ B) (138) 

JC Joo 

Again, for the types of theories usually considered (such as a Klein-Gordon scalar field), the 



surface terms from infinity in eq.(138) vanish. The resulting relation plays an important role in 



defining a natural vacuum state for linear quantum fields in a stationary spacetime [£(]] [21]. 

Consider, now, the case where the nondynamical metric is a flat metric, n ab . We denote the 
(flat) derivative operator associated with 77^ by d a . Let £ a be a translational Killing field of 
so that d a ^ b = 0. Then, clearly, at each point of spacetime the Noether current J associated 
with £ a is linear in the value of £ a at that point. Hence, there exists a unique tensor field, T a b , 
called the canonical energy-momentum tensor, such that 

Ja\...On-\ = ~T e aai...a n _i (139) 

Conservation of J implies conservation of T a b in its first index, i.e., 

d a T\ = (140) 
However, T ab need not be symmetric. Nevertheless, eq.(131) implies that there exists a tensor 



field H abc = H^ c , locally constructed out of 77^ and ip, such that 

jab = jab + Q^cab 

Thus, we have re-derived the well known fact that 7^ always can be "symmetrized" by the 
addition of an identically conserved tensor d c H cab . 

Finally, we note that much of the theory of pseudotensors can be derived by applying the 
results of this Appendix back to the case where the metric again is a dynamical variable in the 



Lagrangian ( 125 ). For a diffeomorphism invariant theory of the type considered in the body of 
this paper, we may introduce a fixed flat metric, rj^, on spacetime, and express the dynamical 
metric, g ab , as, 

9ab = riab + Kb (142) 
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We then may treat n a b and h a b as independent fields, and view our theory as a theory with 
Lagrangian of the form ( |125| ) with the dynamical fields (h a b,ip) and a nondynamical metric 
n a b- One of the (very few) advantages of doing this is that many more quantities qualify 
as "covariant" when rj a b and h a b are viewed as independent fields. In particular, in general 
relativity, no diffeomorphism covariant (n — l)-form, B, satisfying eq,(|80[) can be constructed 
out of g a b, but there is no difficulty in constructing a diffeomorphism covariant B out of the 
independent fields n a b and h a b- Thus, we may change the Lagrangian via 

L -> L' = L - dB (143) 



and still view L' as being of the general form (|125| ) (with r\ a b and h a b viewed as independent 
fields) . Under the change of Lagrangian (|143j) , is modified by 



0' = - ,5B (144) 

(see eq.(|4"2|) above). Consequently, we have t ■ &' = 0, and the canonical energy now is given 
by simply 

£= [ J' (145) 
Jc 



Since the theory has been recast to have a Lagrangian of the form (125) in a spacetime 
with a nondynamical flat metric rj a b, a canonical energy- momentum tensor can be defined by 
eq.( |139| ). We denote this tensor as t a b, and refer to it as a pseudotensor because it depends upon 
the choice of flat metric, r] a b and thus is not covariant with respect to diffeomorphisms which act 
only upon the dynamical fields. For the case of vacuum general relativity with the Lagrangian 
L' of eq.( |143D with an appropriate choice of B, t a b corresponds to the Einstein pseudotensor 
[22]. Note that eq. (p-45|) can be rewritten in terms of t a b as, 



t a b n a t°e 
c 

f j3 



FqcPx (146) 

where the last line holds when C is taken to be the hypersurface t = constant in a global inertial 
coordinate system of r\ a b. 

In order to define a stress-energy tensor corresponding to ( |126| ), we must specify the func- 
tional dependence of the Lagrangian on r] a b for general (non-flat) r] a b- One way to do this would 
be to take L for general rj a b to be given by the substitution (|142| ) in the original Lagrangian. In 
that case, L clearly depends upon n a b and h a b only in the combination ij a b + h a b- Consequently, 
the equations of motion for r\ a b will be satisfied whenever the equations of motion for the dy- 
namical field h a b hold. Note that B need not depend only on the combination r] a b + h a b, so L', 



defined by eq.(143), need not depend only upon this combination. Nevertheless, since addition 
of an exact form to L does not alter the equations of motion, it remains true for L' that the 
equations of motion for rj a b will be satisfied whenever the equations of motion for the dynamical 



field h a b hold. But, this implies that the energy-momentum tensor defined by (126) vanishes by 



virtue of the equations of motion for the dynamical fields. Equation ( |l4l| ) then yields, 

t ab = d c H cab (147) 

This proves that - when the equations of motion are imposed - the pseudotensor t a b always can 
be derived from a "superpotential" H cab . Consequently, the volume integral Ql46j) always can 
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be converted to a surface integral at infinity. This fact, of course, corresponds to our previous 
result, eq.(|83|), which was derived in a much more simple and direct manner. 

When recast in the form ( p.25| ) , the Lagrangian obtained from L' by the simple substitution 
( |142| ) described above will, in general, have a nontrivial, explicit dependence upon the curvature 
of rj a b- However, an alternative procedure for defining a Lagrangian for non-flat 77^ - which 
clearly agrees with L' when rj^ is flat - would be to modify the Lagrangian of the previous 
paragraph by simply setting the terms in L' involving the curvature of 77^ to zero. If we do 
so, the stress-energy tensor defined by ( |126D for this modified Lagrangian will be nonvanishing. 
We denote this stress-energy tensor by t a b- If - as seems plausible - the surface term K does 
not contribute to eq.(p.33|), then the symmetric pseudotensor t a b will be equivalent to t a b insofar 



as the calculation of canonical energy is concerned, i.e., eq.( 146 ) will hold with t a b replaced by 



tab- Note that eqs.( |133| ) and (147) imply that t a ^ also is derivable from a superpotential. Other 
symmetric pseudotensors derivable from a superpotential can be explicitly constructed in the 
case of general relativity (see, e.g., p3| ). 

The dependence of pseudotensors such as t a i, or t a b on the choice of 77^ significantly limits 
their physical interpretation and utility. Indeed, it is difficult to imagine any use to which 
they could be put other than for the definition or calculation of the canonical energy and other 
asymptotic conserved quantities - and this can be accomplished much more straightforwardly 
by the methods described in the body of this paper. Nevertheless, the canonical energy can be 
correctly computed from a pseudotensor via eq.( |l46| ). In particular, if we consider perturbations 
of a stationary solution and choose the flat metric r]^ so that the Killing field, t a , of the 
stationary background is a translational Killing field of r] a b, then 



S 2 £ = J 5 2 t° d 3 x (148) 



is a nontrivial conserved quantity which depends only on the first order perturbation of the 
dynamical fields. For an arbitrary pseudotensor, this formula for 5 2 £ is not very useful because 
to get an expression for the conserved quantity in terms of the first order perturbation, one must 
use the second order field equations to eliminate the terms in 5 2 t°o which involve the second 
order perturbation. However, as shown by Sorkin |22j1 , for the Einstein pseudotensor, the terms 
in 5 2 t°o involving the second order perturbation are separately conserved (irrespective of the 
second order field equations), so the Einstein pseudotensor can be used to obtain a nontrivial 
conserved quantity constructed out of the first order perturbation of the dynamical fields. For 
perturbations of static, electrovac spacetimes in general relativity, this conserved quantity is 
equivalent to the conserved flux integral obtained by Chandrasekhar and Ferrari |24|, |25[| . As 
was shown explicitly in |26| ], the Chandrasekhar- Ferrari conserved current also is equivalent to 
the symplectic current (n — l)-form cv(4>; 5(j), Ct5cj)), which directly yields 5 2 £ by eq. (|136| ) above. 
Note that uj is constructed entirely out of the dynamical fields and their perturbations - in 
particular, no background flat metric need be introduced - so it provides a covariant version 
of the conserved flux integral. However, u is not gauge invariant under infinitesimal gauge 
transformations of the perturbed dynamical fields, so it also does not provide a meaningful 
notion of the local energy density of the perturbation. 
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